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The  production  of  F  and  Cl  atoms  in  an  electrical  discharge  of  F2  or  CI2  has  been  examined  in  a  flow  reactor. 
A  tunable  diode  laser  was  used  to  probe  the  concentration  and  translational  temperature  of  I  atoms  produced 
by  F  and  Cl  atom  reactions  with  HI.  Kinetic  modeling  codes  were  used  to  determine  the  discharge  efficiencies 
from  the  titration  plots  and  the  observed  trends  for  atom  concentration  as  a  function  of  F2  or  CI2  and  pressure. 
These  calculations  indicate  that  the  dc  discharge  used  in  these  experiments  is  100%  efficient  for  F2  flow  rates 
<0.5  mmol  s~*  and  reactor  pressure  <20  torn  The  highest  F2-free  F  atom  flow  rate  that  we  can  generate  is 
1.0  mmol  s~^.  Preliminary  data  for  the  CI2  discharge  indicate  that  this  is  a  much  less  efficient  source  of  Cl  atoms 
with  yields  of  less  than  50%. 


Introduction 

INCH  their  invention  in  the  mid  1960s,  the  most  widely  studied 
class  of  chemical  lasers  has  been  the  hydrogen  halide  systems.' 
These  lasers  operate  on  an  exothermic  energy  release  whereupon 
the  liberated  energy  is  deposited  into  newly  formed  hydrogen  halide 
bonds  to  produce  a  population  inversion  between  vibrational  levels. 
The  following  three-atom  exchange  reactions  typify  the  pumping 
reaction  for  producing  upper  vibrational  levels  in  hydrogen  halides: 


X  •+■  H2 ->  HX(u) -1- H 

(1) 

X-t-HT  HX(u)-f-y 

(2) 

H -1- X2 HX(i)) X 

(3) 

where  X  and  Y  are  the  halogens  F,  Cl,  Br,  I,  and  H  can  be  replaced 
by  deuterium  D.  In  some  cases,  more  than  50%  of  the  available 
exothermicity  is  channeled  into  product  vibration^  and  high  vibra¬ 
tional  levels,  up  to  u  =  1 0,  can  be  populated.  However,  lasing  usually 
occurs  on  A  V  =  —  1  and  —2  transitions  between  the  low-lying  vibra¬ 
tional  levels.  For  example,  the  prominent  rovibrational  laser  lines 
for  the  HF  laser  lie  in  the  (v"  =  2  ->■  u'  =  1)  and  (v"  =  1  ^  u'  =  0) 
bands  between  2.5  and  3.0  fim. 

Recently,  a  population  inversion  and  lasing  were  reported  by 
this  laboratory  on  the  atomic  iodine  transition  at 

1.315  txm  using  the  following  energy  transfer  process^-’: 

NCl(fl'A)  NCIOX^S-)  -t-r^^Pi  j  (4) 

This  system  is  analogous  to  the  well-known  chemical  oxygen  iodine 
laser  (COIL).*  The  NCl(a'  A)  metastable  (To  =  1 . 1  eV)  is  isovalent 
to  singlet  oxygen  in  COIL  and  undergoes  an  efficient  energy  transfer 
reaction  similar  to  the  02(0' A)  -l-I(^P3/2)  reaction  in  COIL.  How¬ 
ever,  in  this  scheme  the  NCl(a'  A)  is  generated  through  a  purely 
gas-phase  chemical  reaction  mechanism 

Cl  -i-  HN3  ->  HCl  N3  (5) 
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Cl-fNs NCl(fl'A)-|-N2  (6) 

and,  as  such,  represents  a  significant  breakthrough  in  chemical  laser 
technology. 

It  should  be  noted  that  reactions  (1-6)  are  preceded  by  a  complex 
set  of  chain  chemistry  where  the  free  radicals  X  and/or  H  are  gen¬ 
erated  in  the  initial  steps.  There  are  many  methods  for  generation 
of  these  radicals,  including  electrical  and  microwave  discharges, 
photolysis,  and  purely  chemical  means  such  as  high-temperature 
combustors.  The  coupling  of  high-enthalpy  flows  found  in  chemical 
lasers  with  these  radical  generator  schemes  produces  large  amounts 
of  heat  in  the  active  medium.  Thus,  extracting  efficient  continu¬ 
ous  wave  laser  operation  involves  managing  homogeneous  reagent 
mixing  times  and  thermal  budgets  of  the  active  medium  with  re¬ 
spect  to  the  kinetics  of  the  system.  Through  optimization  of  several 
flowfield  parameters  such  as  gas  temperature,  velocity,  and  species 
density,  an  appropriate  laser  resonator  and  sub-  or  supersonic  noz¬ 
zle  design  can  be  fashioned  to  accommodate  the  extreme  heat  re¬ 
lease  and  mixing  requirements  found  in  these  systems.  Because 
these  parameters  are  variable,  the  laser  output  power/efficiency  will 
fluctuate  accordingly,  and  a  systematic  choice  of  operational  con¬ 
ditions  is  required  for  efficient  operation.  This  is  especially  true 
of  the  hydrogen  halide  lasers  because  the  vibrational  distribution 
and  small  signal  gain  can  vary  widely  with  operating  conditions. 
The  NCl(a' A)/r(^Pi/2)  laser  system  is  also  sensitive  to  operating 
conditions.  For  example,  to  optimize  the  NCl(fl'A)  and  I*(^Pi/2) 
generation,  complete  dissociation  of  F2  or  CI2  is  essential  because 
IF  and  ICl,  which  are  generated  by 

I  -h  X2  -»■  IX  -I-  X  (7) 

[where  /:f  =  4.5  x  10~'^  and  ka  =  5.5  x  10"'^  cm*  molecule'  s“' 
(Refs.  9  and  10)]  rapidly  quench NCl(a' A)  (Ref.  11)  andr(*Pi/2). 
Note  that  I  (I-l-  CI2)  is  estimated  to  be  the  same  as  k  [I*  (*Pi  /2)  +  CI2] 
(Ref.  10). 

Indeed,  most  chemical  lasers  and  other  high-temperature  react¬ 
ing  flows  contain  complex  chemical  and  fluid  dynamical  processes 
that  require  facUe  temporal  and  spatial  resolution.  Before  the  ad¬ 
vent  of  tunable  diode  lasers,  the  data  required  to  fully  characterize 
these  systems  were  difficult  to  obtain.  However,  a  species’  concen¬ 
tration  and  temperature  can  now  be  readily  characterized  through 
its  spectral  absorption  features  using  commercial  narrow-line  tun¬ 
able  diode  lasers.  The  resulting  spatially  resolved,  high-fidelity,  and 
high-resolution  data  are  particularly  useful  for  characterizing  the 
complex  environment  found  in  chemical  lasers. 

In  this  paper,  we  expand  on  our  earlier  report  of  a  versatile  titra¬ 
tion  technique  that  measures  F  atom  densities  by  using  a  laser- 
based  absorption  technique.'*  The  data  reduction  methodology  has 
been  refined,  and  a  one-dimensional  computational  fluid  mechanics 
analysis  has  been  applied.  Although  this  methodology  is  a  general 
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technique  for  characterizing  F  or  Cl  atom  concentrations  produced 
from  combustion  or  electrical  discharge  devices,  our  specific  inter¬ 
est  is  to  find  conditions  where  large  Fz  and  CI2  free  flows  of  F  or 
Cl  atoms  could  be  generated  using  an  electric  discharge  because 
molecular  halogens  and  interhalogens  are  deleterious  quenchers  of 
NCl(a*A)andr(2Pi/2). 

Experimental  Methods 

A.  +  HI  Atom  Titration 

In  these  experiments,  F  and  Cl  atoms  are  generated  by  a  com¬ 
mercial  electric  (dc)  discharge  (Helios,  Inc.)  in  F2  or  CI2.  Optical 
absoiption  methods  for  determining  absolute  concentration  of  these 
atomic  species  is  difficult  because  the  lowest  energy  allowed  transi¬ 
tions  forF  and  Cl  atoms  lie  in  the  vacuum  The  ^P3/2-^Pi/2 

spin-orbit  transitions'*  for  F  and  Cl,  404  and  881  cm“',  respec¬ 
tively,  can  be  probed  directly,'^''*  but  the  signals  are  quite  weak 
and  a  commercial  source  for  a  404-cm“'  laser  no  longer  exists.  To 
circumvent  these  problems,  titration  methods  have  been  developed, 
where  a  chemical  reagent  that  reacts  with  X  is  added  to  the  flow  to 
produce  a  species  that  can  be  readily  detected  and  converted  into 
absolute  concentration.  For  example,  a  sensitive  technique  for  de¬ 
tecting  F  atoms  in  a  flow  tube  involves  monitoring  the  intensity  of 
HF  (v"  =  3  — >■  u"  =  0)  chemiluminescence  at  880  nm  generated  by 
F  -I-  C2H6  reactions  while  adding  CF3I  as  an  F  atom  titrant.'®  This 
method  is  well  suited  for  small  flows  of  F,  but  would  be  quite  ex¬ 
pensive  and  difficult  to  implement  in  our  large  flow  reactor.  The 
most  common  methods  for  detecting  the  presence  of  Cl  in  a  flow 
are  resonance  fluorescence  and  CI2  afterglow.  Unfortunately,  the 
titrants  CINO  and  C2H3Br  are  not  appropriate  or  cost  effective  for 
our  conditions  (high  flow  rates  and  high  pressure).  In  addition,  these 
methods  do  not  allow  for  simultaneous  concentration  and  tempera¬ 
ture  measurement,  as  well  as  two-dimensional  spatial  resolution  of 
the  F  and  Cl  density  in  the  flow. 

Unlike  most  techniques  that  involve  the  removal  of  a  halogen- 
dependent  chemiluminescent  signal  as  a  function  of  added  titrant 
concentration,  we  monitor  the  generation  of  I  atoms  from  the  re¬ 
action  between  a  halogen  atom  (X  =  F  or  Cl)  and  hydrogen  iodide 
using  an  1.315-/tm  diode  laser  to  probe  the  I(^P3/2)  —  I(^Pi/2)  spin- 
orbit  transition: 

X  -f-  HI  ^  HX  + 1  (8) 

[where  kp  =  4.1  x  10“"  and  ka  =  1.6  x  10"'®  cm^  molecule"'  s"' 
(Ref.  20)].  The  I  atom  absorption  increases  linearly  with  HI  con¬ 
centration  [HI],  until  all  of  the  X  atoms  have  been  consumed.  On 
continued  addition  of  HI,  the  I  atom  absorption  reaches  a  maximum 
and  levels  off  to  a  constant  value.  Ideally,  there  are  no  X  or  I  atom 
loss  processes,  and  the  initial  X  concentration  [X]o  is  given  by  the 
maximum  [I]  observed.  Typically,  however,  at  least  one  or  more 
X  or  I  atom  loss  processes  are  present,  and  [X]o  is  determined  by 
analysis  of  the  I  vs  HI  plots  measured  at  a  variety  of  experimental 
conditions.  In  our  case,  one-dimensional  fluid  mechanics  calcula¬ 
tions  were  used  to  determine  [X]o  and  the  F  and  Cl  atom  production 
efficiency  of  our  dc  discharge. 

B.  Flow  Reactor 

A  schematic  of  the  stainless-steel  flow  reactor  hardware  is  shown 
in  Fig.  la.  The  reactor  is  pumped  by  two  Kinney  850  pumps  in  par¬ 
allel,  each  equipped  with  a  Roots  blower.  The  main  reactor  channel 
is  5.08  cm  wide  by  1.91  cm  high  by  30.5  cm  long  in  the  flow  direc¬ 
tion.  A  2.4-cm-diam  dc  discharge  tube  is  mounted  at  its  entrance.  A 
constant  current  of  286  mA  (~10  kV)  is  applied  to  generate  a  flow 
of  F  or  Cl  atoms  when  a  20%  F2/He  mixture  or  CI2  diluted  with 
He  are  passed  through  the  tube.  The  bulk  He  diluent  is  added  at 
the  interface  between  the  main  reactor  and  the  end  of  the  discharge 
tube.  Hydrogen  iodide  (99%,  Matheson)  diluted  in  He  is  usually 
added  to  the  flow  stream  through  injectors  located  along  the  top 
and  bottom  of  the  reactor.  These  injectors  consist  of  one  row  of 
21  0.081-cm-diam  holes  and  one  row  of  42  0.041-cm-diam  holes. 
The  first  pair  is  located  7.0  cm  downstream  of  the  discharge  exit, 
whereas  the  second  pair  of  identical  injectors  (not  used  in  this  ex¬ 
periment)  is  located  1.8  cm  farther  downstream  and  is  purged  by  a 
small  flow  of  He.  An  axial  ring  injector  located  at  the  exit  of  the 


discharge  tube  is  also  available  for  injection  of  reagents.  The  lo¬ 
cation  of  pressure  transducers  and  thermocouples  is  also  indicated 
in  Fig.  la. 

The  output  of  the  diode  laser  probe  passes  through  Borofloat  win¬ 
dows  (with  a  3-deg  wedge  to  prevent  etalon  effects)  attached  to  both 
sides  of  the  main  reactor  channel.  A  4.25  x  1.91  x  3.05  cm^  vol¬ 
ume  lies  between  the  main  reactor  channel  and  both  inner  window 
surfaces.  To  prevent  the  expansion  of  the  main  reaction  flow  into 
these  areas,  stainless-steel  containment  shrouds  were  installed.  A 
variety  of  shrouds  were  available  depending  on  the  reaction  time  of 
interest.  In  addition,  purge  flows  of  He  were  added  at  the  windows 
and  along  the  shrouds  to  balance  the  main  He  flow.  These  shrouds 
and  purge  flows  not  only  define  the  effective  path  length  for  the 
probe  laser  (5  cm  per  pass),  but  also  protect  the  windows  from  the 
corrosive  chemical  environment  of  the  main  flow. 

The  overall  reactor  pressure  and  linear  flow  velocity  are  controlled 
by  adding  a  large  flow  of  N2  at  the  end  of  the  reactor.  The  linear 
flow  velocity  is  calculated  using 

V  =  MRT/PA  (9) 

where  M  is  the  total  molar  flow  rate  (typically  ~130  mmol  s"'), 
A  is  the  cross-sectional  area  (9.7  cm^),  R  is  the  gas  constant,  P  is 
the  reactor  pressure  (10-30  torr),  and  T  is  the  reactor  temperature, 
determined  from  the  I  atom  probe,  vide  infra.  All  reagent  and  dilu¬ 
ent  flows  were  controlled  by  calibrated  MKS  flow  controllers  or  by 
the  use  of  calibrated  sonic  orifices.  The  upstream  temperatures  are 
monitored  by  type  T  thermocouples  (Omega)  and  the  upstream  pres¬ 
sures  are  controlled  by  a  variety  of  manual  and  electronic  pressure 
regulators  (Proportion  Air). 

C.  Diode  Probe  Laser 

The  optical  path  of  the  diode  probe  is  shown  schematically  in 
Fig.  lb.  The  output  from  a  tunable  1.3145-^tm  diode  laser  (New 
Focus  Lasers  Model  6248;  bandwidth  <500  kHz)  is  optically  sepa¬ 
rated  from  the  rest  of  the  system  by  a  Faraday  isolator.  Beamsplitters 
direct  portions  of  the  output  into  a  Fabry-Perot  etalon  (Tec  Optics 
SA-300)  and  a  heated  iodine  (I2)  cell.  The  I2  cell  temperature  is 
maintained  at  650°C  to  produce  a  known  concentration  of  I  atoms 
that  are  used  to  confirm  the  proper  operation  of  the  laser.  The  main 
portion  of  the  laser  output  is  passed  through  a  half-wave  plate  and 
a  polarizing  beam-splitting  cube.  The  horizontally  polarized  com¬ 
ponent  becomes  the  signal  beam  and  is  passed  on  to  the  reactor, 
while  the  vertically  polarized  component  is  reflected  at  90  deg  and 
is  used  as  the  reference.  The  angle  of  the  half-wave  plate  is  rotated 
to  balance  the  signal  and  reference  beam  intensities.  A  pair  of  In- 
GaAs  detectors  (Physical  Sciences,  Inc.)  monitors  the  reference  and 
signal  intensities.  Before  transmission  through  the  reactor,  the  lin¬ 
early  polarized  signal  beam  is  passed  through  a  quarter-wave  plate 
to  convert  it  to  circular  polarization,  and  a  telescope  focuses  the 
beam  diameter  to  a  1-mm  spot.  A  flat  mirror  loeated  on  the  oppo¬ 
site  side  of  the  reactor  reflects  the  signal  beam  back  on  itself.  The 
return  beam  is  separated  at  the  quarter-wave  plate  and  focused  on 
the  InGaAs  detector.  The  signal  beam  can  be  translated  along  the 
length  of  the  flow  reactor  to  observe  the  temporal  evolution  of  [I] 
and  along  the  perpendicular  axis  to  prohe  the  vertical  profile  of  the 
flow.  The  diode  laser  is  computer  controlled  for  rapid  scans  across 
the  F  =  3  —  F  =  4  hyperflne  component  of  the  I  atom  (5p*  ^Pi/2- 
^Ps/z)  spin-orbit  magnetic  dipole  transition.  Acceptable  signal  to 
noise  levels  are  achieved  by  30-50  scans,  and  I  atom  number  densi¬ 
ties  as  low  as  3  —  5  X  10'^  cm"^  can  be  routinely  detected.  Multiple 
passes  through  the  reactor  can  easily  lower  this  limit:  Densities  as 
low  as  1  X  10®  cm"^  can  he  detected  with  a  2-m  path  length  (pri¬ 
vate  communication,  S.  J.  Davis,  Physical  Sciences,  Inc.,  Andover, 
Massachusetts,  1998). 

The  iodine  atom  number  density  is  derived  from  Beer’s  law: 

I  =  (10) 

The  path  length  L  (two  passes)  is  10  cm,  and  y(v)  is  given  by 

yM  =  {^){x.^Az,,/&jt)fiv)[[n  -  |[/]]  (11) 
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Fig.  la  Schematic  diagram  of  the  flow  reactor;  side  view  of  the  flow  reactor  shows  the  injection  points  for  He  (  rsj  130  mmol  s”*),  Fj/Ch  (0.5- 
1.5  mmol  s“*),  N2,  and  HI  (0.0-6.0  mmol  s~*). 


Hydrogen  Iodide 


Enclosed  Transverse 

Computational  HI  Injector 

Domain  Orifices 


Fig.  Ic  Flow  tube  hardware  and  the  computational  domain  used  in  the  one-dimensional  model;  see  text. 


where  A 34  is  the  Einstein  coefficient  for  the  3,4  hyperfine  transi¬ 
tion  (5.08  s“'),  k  is  wavelength  (1.3145  fitn),  and  f  (v)  is  the  line 
shape  function.  For  our  conditions,  [I*]  =  0,  the  inversion  number 
density  reduces  to  ^  [I],  and  the  ground-state  density  is  obtained  by 
integrating  the  gain  over  frequency; 

'  )/(v)dv  =  — - - ^[I]  (12) 

0  24  Stt 

The  I  atom  spectra  are  fit  by  a  Voigt  function: 
y  =  spectral  area  x  — x  j 


where  wq  and  wl  are  the  Gaussian  and  Lorentzian  linewidths  (full 
width  at  half  maximum),  respectively.  If  collisions  with  He  are  as¬ 
sumed  to  be  the  main  source  of  broadening,  the  I  atom  translational 
temperature  is  calculated  by 

temperature  (AT)  =  M{wq  Ji.\6x  10"^ =  295/u)i  x  3.1  x  Phc 

(14) 

where  vq  is  the  line  center  of  the  transition  (7607.5  cm“’)  and  M 
is  the  atomic  weight  (126.9  g  mol"').  Example  spectra  and  fits  are 
shown  in  Fig.  2. 

Based  on  the  overall  signal-to-noise  ratio  and  the  scatter  of  re¬ 
peated  temperature  measurements  at  the  same  conditions,  the  rel¬ 
ative  error  of  the  temperature  measurement  is  assigned  as  ±30  K. 
The  I  atom  density  is  given  directly  by  the  spectral  area,  and  the 
error  is  quite  small  (<5%). 
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Fig.  2  Sample  I  atom  spectra  showing  F  =  3  —  F  =  4  hyperfine  compo¬ 
nent  of  the  I  atom  (5p®  ^Pi/2-^P3/2)  SP*”  magnetic  dipole  transition; 
signal-to-noise  ratio  is  excellent  and  both  spectra  are  well  fit  by  a  Voigt 
function. 

The  I  atom  flow  rate  is  calculated  by 

i  =  [irVM  (15) 

where  V  is  given  by  Eq.  (9)  and  A  is  the  cross-sectional  reactor 
area.  Note  that  Eq.  (15)  can  only  be  used  if  [I]  is  uniform  along  the 
vertical  axis  of  the  reactor.  Because  the  velocity  depends  on  both 
the  temperature  (±30  K)  and  total  molar  flow  rate,  we  estimate  that 
the  I  atom  flow  rate  is  reliable  to  within  15%.  This  is  consistent 
with  the  systematic  uncertainty  normally  associated  with  flow  tube 
measurements.^' 

D.  Fluid  Mechanics  Simulations:  Methods 

The  model  used  in  this  investigation  to  analyze  the  experiment 
data  is  based  on  the  generalized  one-dimensional  flow  development 
of  Shapiro.^^  This  model  describes  the  flow  through  a  channel  with 
changes  in  area  and  chemical  composition,  secondary  mass  injec¬ 
tion,  wall  friction,  and  heat  exchange  with  the  duct  walls.  The  flow 
changes  are  assumed  to  be  continuous,  and  the  gas  is  assumed  to  be 
real,  with  temperature-dependent  specific  heats.  All  chemical  reac¬ 
tions  are  modeled  with  Affhenius-type  rate  expressions  for  the  for¬ 
ward  rate.  The  equilibrium  constant  for  each  reaction  is  determined 
via  the  Gibb’s  free  energy  method  (see  Ref.  23),  and  the  back¬ 
ward  rate  is  calculated  from  the  forward  rate  and  the  equilibrium 
constant. 

The  primary  aspects  of  the  model  center  on  its  ability  to  predict 
the  changes  in  F  atom  number  densities  by  three-body  recombina¬ 
tion,  wall  recombination,  and  reactions  of  F  and  F2  with  I  and  I2. 
Hence,  these  aspects  of  the  model  are  described  in  detail.  The  finite- 
rate  chemistry  and  reactions  used  by  our  model  are  shown  in  Table 
1  (see  Refs.  24-31).  The  primary  reactions  are  the  production  off 
atoms  via  the  reaction  F  ±  HI,  the  interhalogen  reactions  F  ±  I2  and 
I±F2,  and  the  atom  recombination  reactions.  The  vibrational  dis¬ 
tribution  within  HF  produced  by  the  F  ±  HI  reaction,  the  vibrational 
energy  transfer  reactions  of  HF,  the  energy  transfer  from  HF(i!)  to  I, 


Table  1  F  -t  HI  reaction  rate  constants 


Reaction 

k,  cm^  molecule  '  s  ' 

Reference 

I  atom  production 

F-(-HI->HF(0)-l-I 

1.59  X  10-” 

9 

F-l-HI-»-HF(l)-|-I 

1.91  X  10“" 

— 

F-|-HI->HF(2)-H 

2.39x10-" 

— 

F-fHI^-HF(3)-fI 

2.55  X  10-" 

— 

F-t-HI^-HF(4)-|-I 

3.03  X  10-" 

— 

F-l-HI-^-HF(5)-l-I 

3.98  X  10-" 

— 

F4-HI-^HF(6)-I-I 

2.55  X  10-" 

— 

Total 

1.8  X  lO-'o 

— 

Interhalogen  reactions 

F-|-l2^IF-i-I 

4.3  X  10-'“ 

9 

H-F2^-IF-fF 

1.9  X  lO-''* 

24 

Three-body  recombination 

F  +  F  +  F— >F2  +  F 

1.9  X  10-*T-'exp(629/7’) 

p  _j_  F  4-  F2  F2  +  F2 

1.9x10-*  7-' exp(629/r) 

F-fF+Af  — >■  F2  +  A/ 

7.8  X  10-“exp(629/r) 

H-I-hI^l2  +  I 

5.0x  10-”exp(2120/r) 

9 

I-t-l-l-HI->-l2  +  HI 

3.8  X  io-^“r-“* 

9 

H-I-fl2^l2  +  I 

3.5  X  10-”  exp(2120/r) 

9 

\  +  l  +  M  ^  h  +  M 

5.5  X  10-“‘’exp(575/7’) 

9 

Wall  recombination  reactions 

F  -t  F  -1-  wall  ->■  F2  +  wall 

yi:  =  0.05 

25 

1 4-  H-  wall  ->  I2  +  wall 

yi  =  1.0 

26,  27 

Energy  transfer  reactions 

r-fHF(0)^I-fHF(2) 

1.1  X  10-'2 

28 

r-l-HF(0)-^-H-HF(l) 

1.5  X  10-'* 

— 

r  +  HF(0)^H-HF(0) 

1.8  X  10-'* 

— 

I*-|-F2-^H-F2 

5.0  X  10-'“ 

29 

I*  + 12  I  I2 

3.6x10-" 

30 

P-flF-^-H-IF 

1.3x  10-" 

31 

HFlvl-fHF-^HFlv-l) 

u(5.0x  10-'“T-' 

1 

-1-  HF,  V  =  1-6 

H-5.8xl0-*“r2“) 

HF(v)-l-He^HF(i--  1) 

u(2.6x  io-*“r*) 

— 

-l-He,  u=l-6 

HF(u)  +  F2->HF(v-1) 

0 

1 

0 

X 

— 

-t-F2,  V  =  1-6 

HF(t))4-F-»HF(v-l) 

u[2.7x  10-"exp(-1360/r)] 

— 

-l-F,  V  =  1-6 

HF(t-)-l-HF(u')-»- 

i;'(u±  1)2.7  xl0-"r-' 

— 

HF(u -1- 1)  +  HF(u' -  1), 

V  =  l4,  v'  >v  +  2, 
u'(u+  1)<40 

and  the  deactivation  of  T  and  HF(u)  are  listed  for  completeness,  but 
affect  the  temperature  and  [I]  predictions  by  only  5%.  The  primary 
reaction  F  ±  HI  ->■  HF  ±  I  has  a  nearly  gas  kinetic  rate  constant,  as 
does  the  competing  F  ±  I2  — >•  IF  + 1  reaction.  Proper  accounting  for 
atom  recombination,  particularly  in  the  case  of  I,  is  very  impor¬ 
tant.  The  I  atom  recombination  mechanism  used  here  is  based  on 
the  recommendations  of  Baulch  et  al.®  for  M  =  He,  I2,  and  arbi¬ 
trary  third-body  partners.  The  Arrhenius  rate  constant  expression 
for  three-body  recombination  where  M  =  I  was  extrapolated  from 
the  data  in  Ref.  9  for  both  M  =  I  and  I2.  The  rate  expression  for 
Af  =  HI  was  determined  by  Chang  and  Bums^^  for  the  temperature 
range  200-323  K  and  is  assumed  to  hold  for  the  temperature  range 
of  concern  here  (400-600  K).  For  F  atom  recombination,  we  have 
adopted  the  recommendations  of  the  1977  U.S.  Air  Force  Weapons 
Laboratory  HF/DF  kinetics  review.^^ 

Recombination  of  atoms  at  the  channel  walls  is  modeled  via  a 
source  term  in  the  species  continuity  equation  of  the  form 

Ja  =  yAy/KT\J27TM^[pAw  -  PA2.W  {pA / P/t2)eq] 

where  Ya  is  the  recombination  efficiency  of  species  A,  the  subscript 
W  pertains  to  properties  in  the  gas  at  the  wall,  and  the  subscript  eq 
indicates  equilibrium  conditions  in  the  main  body  of  the  flow.  The 
value  of  Yf  used  here  is  0.05,  the  maximum  value  determined  by 
Kota  et  al.^®  at  room  temperature  for  all  surfaces  investigated  in¬ 
cluding  stainless  steel.  Note  that  the  maximum  value  of  Yf  reported 
by  Kota  et  al.  for  any  surface  was  0.10  at  80  K,  which  indicates 
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that  surface  catalysis  of  F  atom  recombination  is  a  relatively  inef¬ 
ficient  process.  Kota  et  al.  also  report  temperature  dependent  rela¬ 
tions  for  Cl  and  Br  atoms,  (see  Ref.  27).  Tiie  I  atom  recombination 
coefficient  used  in  this  model  is  1.0,  based  on  the  recommenda¬ 
tions  of  Perram  and  Hager^®  and  the  work  of  Fisk  and  Hays.^ 
This  value  is  also  consistent  with  the  recombination  coefficients 
recently  reported  by  Kota  et  al.  for  Cl  and  Br  atoms^"*  (also  see 
Ref.  26). 

The  computational  domain  for  the  physical  space  within  the  ex¬ 
periment  simulated  with  this  model  extends  from  the  exit  of  the 
discharge  tube  to  the  point  at  which  the  I  atom  measurements  ate 
taken  in  the  streamwise  direction,  from  wall  to  wall  and  shroud  to 
shroud  in  the  lateral  direction,  and  from  wall  to  wall  in  the  vertical 
direction,  see  Fig.  Ic.  The  discharge  tube  is  treated  as  a  black  box 
entity,  for  which  only  the  exiting  mass  flow  is  known;  the  pressure, 
temperature,  and  the  F2  dissociation  fraction  of  the  flow  exiting 
the  discharge  tube  are  unknown  quantities  that  are  determined  by 
matching  the  model  predictions  to  measured  quantities  downstream. 
For  the  range  of  pressures,  temperatures,  and  mass  flow  rates  con¬ 
sidered  here,  the  flow  exiting  the  discharge  tube  is  subsonic,  and 
in  all  simulations  of  the  experiment,  the  flow  is  treated  as  subsonic 
throughout. 

An  important  simplifying  assumption  of  our  calculations  is  that 
the  Hi/He  flow  is  modeled  with  instantaneous  mixing  at  the  same 
streamwise  position  relative  to  the  discharge  tube  exit  at  which  the 
actual  flow  is  injected.  Secondary  He  is  injected  at  the  position  of 
the  second  row  of  injectors,  consistent  with  the  experiment.  Because 
the  pressure  across  the  injector  orifices  is  always  greater  than  2:1, 
the  sonic  conditions  are  used  to  specify  the  temperature  and  pressure 
for  each  injectant  flow. 

The  numerical  technique  used  to  solve  the  governing  equations 
of  the  model  is  a  first-order  Adams-Bashfourth-Moulton  method 
(see  Ref.  35).  In  this  numerical  integration,  the  nominal  step  size 
Ax  is  0.05  cm.  This  value  was  chosen  by  iteration  until  changes  in 


a)  laser  probe  position 


b)  laser  probe  position 

Fig.  3  Vertical  profiles  at  various  positions  along  the  reactor;  I  atom 
laser  probe  can  be  translated  along  the  transverse  and  vertical  axis  of 
the  reactor.  Profiles  are  flat  for  d  >  5.S  cm. 


the  integrated  solution  were  <  1  %.  The  step  size  after  gas  injection 
locations  was  decreased  to  Ax/200  and  returned  to  Ax  over  a  dis¬ 
tance  of  4 Ax.  The  solution  is  marched  through  the  channel  in  the 
flow  direction  giving  only  upstream  dependence  and  not  directly 
capturing  the  elliptic  nature  of  the  subsonic  flow.  The  calculated 
downstream  pressure  is  determined  by  varying  the  inflow  pressure 
until  the  model  prediction  and  measurement  agree.  This  technique 
is  also  used  to  match  the  experimentally  determined  temperatures. 
Reactant  flow  rates  are  input  to  the  model  as  fixed  quantities  and  are 
not  affected  by  the  pressure/temperature  iterations. 

Each  titration  represents  a  series  of  experiments  where  all  param¬ 
eters  and  flow  rates  are  held  constant  except  [HI] .  As  the  HI  flow  rate 
approaches  0,  the  heat  release  of  the  F  -|-  HI  reaction  is  minimized, 
and  the  flow  conditions  approach  those  for  the  flow  without  HI. 
Hence,  the  initial  temperature  exiting  the  discharge  is  determined 
by  matching  the  temperature  and  pressure  with  the  model  for  the 
lowest  HI  flow  rate  in  the  titration.  On  the  other  hand,  because  the 
F2  dissociation  fraction  exiting  the  discharge  tube  is  also  unknown, 
this  quantity  is  determined  by  matehing  the  I  atom  number  density 
measured  for  the  highest  HI  flow  in  each  titration  series.  The  high 
[HI]  data  are  used  because  they  have  the  highest  signal-to-noise  ra¬ 
tio  and  are  the  least  kinetic  model  dependent.  With  the  dissociation 
fraction,  pressure,  and  temperature  fixed,  the  model  is  run  for  the 
remaining  HI  flow  rates  in  the  titration. 


a)  Solid  and  broken  lines  indicate  the  model  calculation  results 


Fig.  4  Typical  F  +  HI  titration  plots,  I  atom  densities,  I  atom  flow  rates, 
and  temperatures  as  a  function  of  added  HI  at  F2  =  0.45  mmol  s~^: 
•,  P  =  12  torr;  ■,  P  =  15  torr;  and  A,  P  =  20  torn 
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a)  Solid  and  broken  lines  indicate  the  model  calculation  results 


b)  Solid  line  indicates  tbe  theoretical  limit  for  100%  dissociation 


Fig.  5  Typical  F  +  HI  titration  plots,  I  atom  densities,  I  atom  flow  rates, 
and  temperatures  as  a  function  of  added  HI  at  F2  =  1.30  mmol  s~^: 
•,P=  12  torr;  ■, /*=  15  torr;  and  A,f’=  20  torn 


Experimental  Results 

A.  Vertical  Profiles 

An  important  advantage  of  our  experimental  technique  is  the  abil¬ 
ity  to  measure  I  atom  concentrations  and  temperatures  along  both 
the  transverse  and  streamwise  directions  of  the  flow  reactor.  Fig¬ 
ure  3a  shows  the  vertical  profiles  of  the  I  atom  flow  rate  measured 
at  F2  =  0.45  mmol  s“'  with  excess  HI,  F  =  15  torr,  and  five  dif¬ 
ferent  positions  along  the  reactor.  A  uniform  concentration  profile 
should  be  exhibited,  indicative  of  a  constant  velocity  distribution 
along  the  vertical  axis.  The  data  measured  2.5  cm  from  the  HI  inlet 
are  peaked  in  the  center  and  fall  off  toward  the  edges.  Clearly,  a 
uniform  flowfield  has  not  been  established  and  these  data  are  not 
useful  for  determining  [I]  and  T.  The  presence  of  at  least  one  I 
atom  loss  process  seems  evident  by  the  factor  of  2  loss  of  I  atoms 
between  d  —  5.5  and  15.5  cm.  However,  subsequent  chemilumi¬ 
nescent  measurements  indicated  significant  flaring  of  the  flow  and 
leakage  around  the  10.5-  and  15.5-cm  flow  containment  shrouds. 
Consequently,  only  the  5.5  cm  datum  is  considered  to  be  reliable. 
Future  reactor  designs  will  resolve  this  problem.  Figure  3b  shows  the 
I  atom  temperature  profiles.  Ald  =  2.5  cm,  the  temperature  peaks 
in  the  center  and  is  ~100  K  lower  at  the  edges.  The  I  atom  flow 
data  violate  mass  balance,  which  indicates  that  the  plug  flow  has 
not  been  established.  As  the  distance  between  the  reagent  inlet  and 
the  observation  point  increases,  the  temperature  decreases  and  the 
temperature  profiles  are  no  longer  peaked. 


B.  F  +  HI  Titrations 

The  data  from  six  typical  titrations  are  shown  in  Figs.  4  (0.45 
mmol  s"'  F2)  and  5  (1.3  mmol  s"'  F2).  Figures  4a  and  5a  show  [I] 
plotted  vs  the  flow  rate  of  added  HI  at  three  different  pressures  and 
d  =  5.5  cm.  In  all  cases,  the  titration  plots  have  both  a  linear  and  a 
plateau  region.  Solid  and  broken  lines  indicate  the  model  calculation 
results  for  [I].  Figures  4b  and  5b  show  the  corresponding  I  atom 
flow  rates  as  a  function  of  added  HI,  relative  to  the  theoretical  limit 
(indicated  by  the  solid  line),  which  assumes  100%  dissociation  of  F2, 
1;1  stoichiometry,  and  no  F  or  I  atom  losses.  Finally,  the  measured 
temperatures  are  plotted  Figs.  4c  and  5c  and  are  compared  to  the 
average  values  from  the  model  calculations. 

In  all  but  one  case,  the  model  accurately  reproduces  the  exper¬ 
imentally  measured  [I];  the  only  exception  is  for  F  =  12  torr  and 
F2  =  0.45  mmol  s"’,  where  the  model  underestimates  [I].  However, 
as  Fig.  4b  shows,  the  I  atom  flow  rate  at  12  torr  also  exceeds  the 
theoretical  limit.  Although  this  may  be  indicative  of  the  flowfield 
problem  described  earlier,  the  error  is  only  ~20%,  and  the  problem 
is  not  serious.  In  general,  however,  the  agreement  between  the  model 
and  the  experiment  is  very  good,  with  the  values  and  slopes  of  the 
titration  curves  matched  very  well  in  most  cases.  The  agreement  be¬ 
tween  the  model  and  experiment  for  the  lower  HI  flow  rates  of  each 
titration  suggests  that  the  mechanisms  that  interfere  with  a  pure  titra¬ 
tion  of  F  atoms  (i.e.,  I  -f  F2,  F-l- 12,  three-body  recombination,  and 
wall  recombination)  are  modeled  accurately.  The  average  calculated 
temperature  at  5.5  cm  and  F2  =  0.45  mmol  s"*  is  418  K,  consistent 
with  the  experimental  results.  The  data  for  F2  =  0.45  mmol  s“'  are 
also  consistent  with  100%  dissociation.  The  dissociation  efficiency 
is  less  than  70%  for  the  data  with  F2  =  1.30  mmol  s“*  (see  Fig.  5). 
The  calculated  temperature  for  F2  =  1 . 3  mmol  s“  *  is  higher  than  for 
F2  =  0.45  mmol  s“’ ,  Tav  ^aic  =  488  K,  and  is  consistent  with  the  mea¬ 
sured  values.  The  agreement  between  the  predicted  temperatures  via 
the  line  shape  analysis  and  the  model  predictions  tends  to  validate 
both  the  line  shape  analysis  and  the  model. 

Figures  6  and  7  summarize  the  titration  results  as  a  function  of  F2 
flow  rate  and  pressure.  In  all  three  panels,  the  data  points  represent 
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Fig.  7  Dissociation  fraction  vs  pressure;  Fi  dissociation  fraction,  a  = 
Iitiax/2*F2,  is  plotted  as  a  function  of  the  reactor  pressure. 


Fig.  8  Typical  Cl  +  HI  titration  plots;  I  atom  flow  rates  and  tempera¬ 
tures  as  a  function  of  added  HI  at  a  variety  of  experimental  conditions. 


the  average  I  atom  flow  rate  from  the  plateau  region  of  the  corre¬ 
sponding  titration  curves.  Figure  6  demonstrates  the  effect  of  in¬ 
creased  F2,  and  although  the  I  atom  flow  rate  is  shown  to  increase 
with  larger  F2  flow  rates,  the  I  atom  yield  given  by 


a  =  I  atom  yield  = 


I  atom  flow  rate 
2*F2  flow  rate 


(17) 


is  inversely  proportional  to  the  F2  flow  rate.  These  data  indicate  that 
both  the  F2  flow  rate  and  reactor  pressure  affect  the  performance  of 


the  dc  discharge.  The  highest  F2-free  flow  rate  of  F  atoms  that  we 
are  capable  of  generating  is  1.0  mmol  s"'.  As  shown  in  Fig.  7,  for 
F2  <  0.44  mmol  s” '  and  P  <20  torr,  1 00%  dissociation  is  achieved. 
At  higher  F2  flow  rates  and  pressures,  the  dissociation  efficiency 
decreases  linearly. 

C.  Cl  +  HI  Titrations 

A  set  of  typical  Cl  -f  HI  titrations  are  shown  in  Fig.  8.  The  atom 
flow  rate  generated  by  the  CI2  discharge  is  significantly  smaller  than 
that  obtained  from  F2.  In  addition,  the  temperature  is  slightly  lower, 
T  ~  400  K.  Increasing  the  CI2  flow  rate  does  not  generate  more 
Cl  atoms  and  increasing  the  pressure  also  degrades  the  discharge 
performance.  The  highest  Cl  atom  flow  rate  that  we  can  generate 
is  0.3  mmol  s“'  independent  of  the  initial  flow  rate  of  CI2,  and  the 
average  temperature  is  ~380  K.  The  best  dissociation  efficiency, 
a  =  50%,  is  achieved  for  CI2  =  0.5  mmol  s'"'  and  P  =  15  torr.  This 
cursory  analysis  clearly  indicates  that  CI2  is  inferior  to  F  -|-  HCl/DCl 
as  a  Cl  atom  source.  As  a  result,  no  model  calculations  were  per¬ 
formed. 


Conclusions 

The  present  results  demonstrate  the  utility  of  the  diode  laser  as  a 
probe  of  atom  flow  rates  and  temperatures.  Based  on  the  experimen¬ 
tal  data  and  one-dimensional  fluid  mechanical  modeling  results,  we 
have  established  the  efficiency  of  our  dc  discharge  for  generation  of 
laige  F  and  Cl  atom  densities.  The  X  -H  HI  titration  reactions  effec¬ 
tively  convert  the  F  and  Cl  atoms  to  iodine  atoms,  which  are  detected 
by  near  IR  laser  absorption.  This  method  has  clear  advantages  over 
alternative  techniques  such  as  resonance  fluorescence  or  absorption, 
CI2  recombination  afterglow,  and  HF  ( An  =  3)  chemiluminescence. 
Specifically,  the  laser  probe  technique  gives  both  atom  densities  and 
temperatures  and  can  give  two-dimensional  spatial  resolution  when 
the  beam  is  translated  along  the  vertical  and  streamwise  axes  of  the 
reactor.  When  the  path  length  is  increased,  extremely  low  concen¬ 
trations  can  be  measured. 
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